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The  experimental  and  theoretical  work  performed  under  this  Contract  was 
concerned  with  a scale  model  study  of  various  properties  of  Arctic  ice,  in 
particular  the  reflectivity  of  sonic  signals  from  the  sea-water/ice-cover 
interface.  The  scaling  was  accomplished  hy  considering  the  usual  thickness 
of  the  Arctic  ice  cover  and  the  wavelength  of  the  usual  sonar  signal,  and  then 
changing  the  ice  thickness  to  a few  millimeter  and  the  frequency  of  the  sonic 
signal  to  a few  megahertz;  this  scaling  maintains  the  same  value  of  the  product 
(frequency  times  thickness)  for  both  cases,  the  situation  in  the  Arctic  and 
the  scaled  version  of  the  phenomena  investigated  in  the  laboratory. 

Although  reflectivity  coefficients  can  be  easily  calculated  from  time- 
honored  formulas,  these  calculations  do  not  consider  the  fact  that  plane  waves 
of  infinite  extent  are  not  possible  in  nature,  and  that  all  signals  do  have  a 
finite  dimension,  i.e.,  one  invariably  deals  with  a finite,  bounded  beam  of 
sonic  signals.  This  departure  from  the  simnlistic  view  creates  new  problems  in 
the  evaluation  of  reflectivity  , and  thus  a great  part  of  the  work  performed  was 
concerned  with  the  formulation  of  appropriate  expressions  describing  the  actual 
situation.  This  formulation  vas  then  experimentally  verified  for  easily  handled 
interfaces,  namely  nonmelting  plates  immersed  and  also  floatinr  on  water.  This 
was  done  initially,  before  experiments  with  ice  plates  were  conducted,  in  order 
to  check  out  the  basic  premises  of  the  newly  developed  theoretical  approach  to 
the  problem.  The  next  step  was  to  check  into  the  angular  dependence  of  the  in- 
cident/reflected beam  ratio,  and  it  was  found  that  here  too  the  plane  wave  theory 
leads  to  useless  results,  and  the  bounded  beam  approach  had  to  be  investigated. 
This  was  again  accomplished,  first  for  ice-simulating  plates  and  water  interfaces 
and  then  for  ice  plates  of  millimeter  thickness. 

In  the  course  of  these  investigations  it  was  realized  that  the  initial  am- 
plitude of  the  impinging  signal  enters  the  description  of  the  reflection  if  non- 
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linear  processes  occur  during  the  proportion  in  the  water  or  upon  re- 
flection. These  frequency  dependent  and  amplitude  dependent  phenomena  became 
increasingly  important  and  non-negligible  in  the  course  of  the  work  that  these 
processes  warranted  a separate  investigation  by  themselves.  Thus , with  the 
agreement  of  the  Code  sponsoring  the  present  Contract,  a renewal  proposal  was 
submitted  to  Code  1*21  for  the  expressed  purpose  to  investigate  the  parameters 
of  nonlinearity,  in  a most  general  manner,  before  a continuation  and  conclusion 
of  the  present  contract  could  be  considered  in  the  previously  anticipate  manner. 

nevertheless,  the  present  Contract  has  produced  a number  of  results  we  had 
anticipated,  and  these  results  are  listed  in  this  Report.  The  present  Contract 
actually  consisted  of  two  seperate  activities:  Contract  IJ000ll*-69-A-O22  0-007 , 
sponsored  through  Code  1*15,  and  upon  reorganization  within  OUR,  Contract 
N0001U-75-C-0333,  sponsored  through  Code  l*6l.  Although  this  reorganization  and 
numbering  of  the  Contracts  did  not  imply  different  efforts,  it  caused  a Final 
Report  to  be  submitted  at  the  time  of  the  reorganization.  Clearly,  the  work 
was  not  completed  in  all  of  its  phases  at  that  time,  thus  the  first  Final  Report 
described  the  work  accomplished  until  that  time.  Various  Technical  Reports 
have  been  submitted  throughout  the  periods  covered  by  both  Contracts,  and  thus 
this  Report  will  refer  to  prior  submissions  by  means  of  citing  the  titles  and 
abstracts  of  the  various  publications  originating  from  this  activity  in  the 
Contract  period. 

The  more  recent  results,  not  covered  in  previous  Ttc..:.ical  Reports  will  be 
discussed  in  more  detail.  In  this  connection  it  should  be  noted  that  an  exten- 
sion of  the  present  Contract  {no  cost)  was  requested  and  this  request  had  been 
granted  because  the  final  reduction  of  the  rather  complex  end  involved  computer 
results  generated  and  the  comparison  of  experimental  results  obtained  in  the  lab 
and  presently  stored  or.  Video  Tepe  have  been  most  time-consuming  end  have  necessi- 
tated computation  adjustments  previously  thought  to  be  minor.  Thus  the  present 


report  is  a Final  Report  to  which  an  Addendum  will  be  submitted  as  soon  as  the 


data  reduction  and  comparison  is  completed. 

Section  I of  this  Report  contains  copies  of  titles  and  abstracts  of  work 
performed  and  already  submitted  to  the  sponsoring  Code.  Some  comments  concerning 
continuity  and/or  reasons  for  performing  the  listed  vork  are  included  in  this 


section.  Section  II  describes  the  work  not  reported  on  in  Section  I and  accom- 
plished up  to  the  expiration  of  the  Contract. 


i ..  • 


Section  I 


The  preliminary  findings  during  the  first  few  months  of  the  Contract  were 
concerned  with  an  overall  investigation  of  the  parameters  which  play  a role  in 
the  determination  of  ice  thickness  via  reflection  from  the  Arctic  ice  cover. 
These  findings  were  reported  on  in  Technical  Report  IJo.  1. 


SCALE  MODEL  ULTRASONIC  STUDY  OF  ARCTIC  ICE 
TECHNICAL  REPORT  NO.  1 


Backscattering  of  ultrasonic  signals  from  solid  plates  is  interpreted 
to  yield  information  about  the  flexural  modes  of  the  plate  and  the  son- 
ic velocities  in  sea  ice  as  well  as  the  fractional  ice  content  of  the 
water/sea-ice  layer. 


It  became  apparent  that  the  shear  velocity  in  ice  influences  the  reflection 
properties  of  the  ice  to  a great  extent,  and  since  this  velocity  is  ordinarily 
of  no  great  interest  in  other  Arctic  work,  little  information  was  available  about 
it.  Therefore,  as  a pre-requisite  for  completing  the  work,  an  investigation  was 
undertaken  which  resulted  in  the  following  paper: 


Determination  of  Transverse  Wave  Velocity  in  Transition  Layer 
of  Sea  Ice  From  Reflection  of  Water-Borne  Sound 

Measurement  of  reflection  coefficients  at  water-sea  ice  in- 
terfaces were  reported  by  Langleben  11970],  and  the  values  of 
the  sound  velocity  in  the  skeleton  layer  as  well  as  its  density 
were  deduced  from  the  data.  Langleben  used  a formula  given 
by  Rayleigh  [1945]  that  is  applicable  when  the  interface  is 
formed  by  two  liquids.  The  use  of  this  formula  resulted  in  the 
finding  that  the  density  of  the  skeleton  layer  ice  was  com- 
paratively high  (0.996  g cm'*)  and  the  sound  velocity  was  low 
(1810  m/s).  It  was  also  noted  by  Langleben  that  for  large 
angles  of  sound  incidence  the  reflection  coefficients  were  con- 
siderably lower  than  the  value  of  unity  predicted  by  the 
Rayleigh  formula.  These  data  were  thus  omitted  from  the 
evaluation. 

A reevalualion  of  the  data  is  presented  here  that  removes 
some  of  the  difficulties  encountered. 


It  was  found  advantageous  to  have  available  a method  with  which  the  sonic 
velocities  in  ice  or  any  other  solid  can  be  measured  with  one  technique,  since 
both  of  these  velocities  are  required  as  data  in  the  evaluation.  This  method 
was  developed  and  was  described  in: 


A Simple  Technique  for  Simultaneous  Measurement 
of  Longitudinal  and  Shear  Wave  Velocities  of  Solids 


Sound  velocities  in  solids  can  be  determined  indirectly  by  In  this  paper  we  | ropoae  a simple  method  in  which  the 

measuring  ultrasonic  reflection  or  transmission  at  liquid-  existence  of  critical  angles  of  the  solid  is  not  a requirement 

solid  boundaries.  The  incidence  of  the  ultrasonic  beam  may  nor  is  it  neerssary  to  use  two  types  of  transducers.  With 

be  normal  to  the  surface  of  the  bulk  solid  in  which  case  this  method,  both  the  longitudinal  and  the  shear  wave 

only  the  longitudinal  sound  velocity  can  be  determined  velocities  in  the  solid  are  measured  simultaneously  via 

from  measurements  of  the  reflection  coefficient,  or  the  Schlieren  techniques  using  the  same  incident  longitudinal 

sound  beam  may  impinge  from  the  liquid  at  oblique  in-  sound  beam, 

cidence  in  which  case  one  can  obtain  the  sound  velocities  in 
the  solid  hy  observing  angles  of  maximum  reflection. 


In  connection  with  this  investigation  it  was  soon  found  that  reflectivity 
at  the  critical  angles  of  incidence  (primarily  the  Rayleigh  angle)  is  markedly 
different  froc  the  reflection  at  other,  noncritical  angles.  A study  of  these 
drastic  changes  upon  reflection  can  best  be  done  with  polycrystalline  solids 
(not  restricted  to  ice) , and  this  stucty  resulted  in  a paper  listed  below: 


second  Technical  Report: 


Sonic  Reflectivity  from  Sea-Ice/Water  Interface 
Technical  Report  No.  2. 


Reflectivity  curves  are  calculated  for  sea-ice/water  boundaries  in 
which  the  densities  and  sonic  velocities  in  the  two  media  are  chan- 
ged in  increments.  The  resulting  changes  in  the  reflectivity  curves 
are  interpreted  in  terms  of  the  influence  of  the  parameters.  Anal- 
ysis is  given  indicating  under  which  conditions  reflectivity  measur- 
ements can  or  cannot  yield  information  about  the  values  of  sonic  vel- 
ocities and  densities  of  the  sea-ice  and  the  water.  Possible  mis- 
leading measurements  are  indicated  and  the  influence  of  Rayleigh 
waves  on  measurement  is  described. 


At  this  point  it  beceme  apparent  that  nonlinear  effects  in  the  water  or  the 
ice  will  influence  the  results  if  the  frequency  and/or  the  amplitude  of  the  pro- 
bing sonic  signal  are  high.  In  order  to  be  eble  to  eliminate  this  unwanted  effect 
it  became  necessary  to  understand  the  parameters  which  are  responsible  for  the 
presence  of  the  effect.  Although  the  effect  initially  was  to  be  excluded  from 
the  considerations  its  mere  presence  made  an  investigation  of  it  mandatory  before 
its  elimination  could  be  guaranteed. 

Since  the  best  way  of  studying  the  effect  is  contained  in  an  examination 
of  surface  waves  on  single  crystals  (due  to  its  orientation  dependent  values  of 
the  surface  wave  velocities,  unfortunately  absent  in  polycrystalline  substances 
like  large  quantities  of  ice),  we  conducted  further  experiments  on  nonlinear 
surface  wave  interactions.  The  results  were  published  in 


Further  investigation  of  noncollinear  surface-wave  interactions 


Two  noncollinear  surface  waves  of  5 and  7 MHz  have  been  exciled  on  the  surface  of  single-crystal  copper 
samples  in  such  directions  that  the  waves  could  interact  via  the  nonlineanties  of  the  medium  to  produce  a 
third  surface  wave  at  the  sum  frequency  of  12  MHz.  Using  optical  techniques,  the  three  waves  involved  have 
been  detected  on  two  crystal  surfaces,  the  (001)  and  (III)  planes  of  copper,  and  the  amplitude  of  each  of  the 
waves  has  been  measured. 


ry- 

I These  results  and  other  obtained  outside  the  work  related  to  the  present 

Contract  formed  the  basis  for  the  renewal  proposal  alluded  to  above. 

Clearly,  any  deep  re-investigation  of  the  nonlinear  effects  would  have 
brought  the  efforts  outside  the  task  of  the  Contract,  thus  we  concentrated  on 
the  continuation  of  reflectivity  work,  taken  into  account  all  precautions 
necessary  to  avoid  nonlinear  behavior  of  the  substances  involved. 

It  was  found  that  the  loading  of  the  plate  does  change  the  value  of  the 
surface  wave  velocity,  and  this  velocity  is  an  important  quantity  in  the  cal- 
culation of  reflectivity.  Ultimately,  the  surface  wave  velocity  can  be  deter- 
mined from  an  observation  of  the  reflectivity  behavior,  procipally,  the  phase 
of  the  reflection  at  incidence  at  the  Rayleigh  angle.  The  work  on  the  phase  of  • 
the  reflection  was  described  in: 

Phase  of  ultrasonic  reflection  at  Rayleigh  angle  incidence 

The  phaic  of  ui  ultrasonic  bam  reflected  from  • liquid-solid  interface  at  and  nut  Rayleigh  angle 
incidence  U measured  using  a new  technique  Results  indicate  that  two  waves  comprise  the  reflected 
beam,  te.,  a specularly  reflected  wave  and  a rcradiated  surface  wave,  which  do  nol  change  phase  as 
the  angle  of  incidence  variea,  which  maintain  a 180*  phase  difference  hetv-een  them,  and  which 
propagate  collincarly  only  at  the  Rayleigh  angle. 


Approximating  or  simulating  the  properties  of  ice  as  closely  as  possible 
by  using  plastic  plates  in  the  inital  experiments  on  the  reflectivity  showed 
that 

a)  the  finite  thickness  of  the  plate  and  the  frequency  of  the  sonic 
sif^ial  used  results  in  the  existence  of  more  than  one  critical  angle 
of  reflection  where  reflectivity  data  behave  anomalously, 

b)  the  reflection  coefficient  for  bulk  substances  is  practically  unrelated 
to  reflectivity  from  liquid/solid/liquid  loyer  geometries. 

These  points  were  discussed  at  the  London  IGA  under  the  title: 


””  1 
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EIGHTH  INTERNATIONAL  CONGRESS 
ON  ACOUSTICS.  LONDON  1974 

ULTRASONIC  RERADIATION  FROM  RAYLEIGH  AND  LAMB  WAVES 
AT  LIQUID-SOLID  INTERFACES 

Moyer  W G Georgetown  University,  Washington  ,D  C /USA 

Plena  T J 
Behravesh  M 
INTKO DUCT ION 

The  surface  wave  excited  by  a sound  beam  Incident  on  a liquid-solid 
bulk  Interface  at  the  Rayleigh  angle  reradlates  energy  Into  the  liquid 
at  the  Rayleigh  angle.  The  reradlatlon  and  specular  reflection  prod- 
uce a resultant  sound  field  partially  described  by  Schoch's  "beam  dis- 
placement" (1),  and  recently  more  completely  by  Bertonl  and  Ifcmlr  (2). 

In  a similar  manner,  Iamb  waves  can  be  excited  at  liquid-solid 
plate  Interfaces  to  produce  similar  "beam  displacements"  (3). 

THEORETICAL  CONSIDERATIONS  AND  EXPERIMENTAL  RESULTS 

Phase  velocities  V of  free-plate  Iamb  modes  depend  on  the  product  fd 
(frequency  x plate  thickness)  and  the  longitudinal  and  shear  veloci- 
ties of  the  solid.  The  lowest  l'ew  I/u-b  modes  are  shown  for  brans  in 
Pig.  1,  and  for  Plexiglas  in  Fig.  2. 

If  the  plate  is  immersed  in  water,  loading  changes  the  V.  New  val- 
ues are  found  by  adding,  to  the  free-plate  equations,  terms  containing 
the  ratio  r « [density  of  water/density  of  the  solid].  V can  be  ex- 
pressed in  terras  of  the  angle  of  sound  incidence  a,  from  Snell's  law. 

'•  Calculations  for  brass  (p  » 8.6)  show  that  the  values  of  V differ 
very  little  for  free  or  for  loaded  plates  (Fig.  1),  while  for  solids 
with  p « 1,  the  two  sets  of  curves  differ  markedly  (Flfi:.  i and  3)- 

Schlicr-o  techniques  are  used  to  experimentally  verify  the  results. 
For  a given  solid  and  a fixed  fd,  those  a were  determined  where  "beam 
dlsplaceracnt"occurred,  giving  values  of  V.  Results  (Figs.  2 and  3) 
show  that  a free-plate  approach  should  not  be  used  for  plastic,  ice, 
and  ether  low-density  plates  in  water. 


Acknowledgment:  fupported  by  the  Office  of  Naval  Research,  U.S.  Navy. 
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IiOO0ll;-69-A-022O-0O7.  In  terms  of  the  tesk,  there  was  no  discontinuity  in  effort 
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be  subiuitted  at  the  conclusion  of  the  computation  evaluation.  This  evaluation  is 
now  being  accomplished  at  no  additional  cost  to  the  Contract. 
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Section  II 


As  was  pointed  out  above,  the  determination  of  ice  parameters,  like 
thickness,  via  sonic  reflectivity  depends  on  the  appropriate  combination  of 
information  about  the  reflection  properties,  in  general,  of  a bounded  beam 
from  a solid/liquid  interface,  and  on  the  reflectivity  parameters  of  ice  in 
its  various  mechanical  forms  (i.e.,  density,  sonic  velocities).  Thus  at  the 
start  of  the  second  Contract,  the  effort  was  two-fold:  determine  standard 
reflectivity  for  plane  waves  from  an  infinite  half  space  combination  ice/sea 
water,  and  to  continue  the  investigation  of  bounded  beam  reflectivity  in  gene- 
ral. The  first  topic  was  investigated  and  the  results  are  given  in  the  follow- 
ing: 

Determination  of  sonic  velocities  from  reflectivity  losses  at  sea 
ice/water  boundaries 


Refleclivhies  ire  calculated  as  functions  of  sonic  incidence  angle,  densities,  and  sonic  velocities.  The 
resulting  sets  of  curves  are  used  to  determine  under  which  conditions  reflectivity  measurements  can 
or  cannot  yield  information  about  the  values  of  sonic  velocities  and  density  of  sea  icc. 


This  paper  describes  the  various  differences  one  may  expect  whenever  a sonic 
signal  is  reflected  from  the  ice/sea  water  boundary.  It  does  yield  information 
about  the  angular  ranges  where  the  three  important  parameters  (density,  longitu- 
dinal, shear  velocity)  influence  markedly  the  plane  wave  reflection  coefficient, 
and  where  variations  in  any  of  these  parameters  do  not  significantly  change  the 
reflection  coefficient  magnitude. 

The  second  topic  which  was  investigated  does  relate  to  the  problem  of  finite 
beam  reflection  coefficient  versus  plane  wave  reflection  coefficient,  and  as  a 
subdivision  of  these  investigations,  the  influence  of  solid  geometry  on  the  finite 
beam  width  reflection.  The  geanetry  in  turn  can  be  divided  into  cases  where 
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the  plate  reflector  is  bounded  by  the  seme  liquid  on  both  sides,  by  dissimilar 
liquids  on  the  two  sides  of  the  solid  plate,  and  by  a liquid  on  one  side  and 
vacuum  or  eir  on  the  other. 

The  first  investigation  in  this  regard  was  concerned  with  a comparison  of 
solid  bulk  versus  solid  plate  reflectors  and  the  types  of  surface  waves  which  can 
be  set  up  in  the  reflection  process.  A description  of  these  surface  waves, 
Rayleigh  or  Lamb  waves  is  contained  in 


Rayleigh  and  Lamb  waves  at 
liquid— solid  boundaries 


The  equations  governing  Rayleigh  and  Lamb  mode  propagation  are  examined  for  free 
and  for  liquid-loaded  solids.  Examples  are  given  to  show  under  what  conditions  the 
free-solid  approach  yields  acceptable  solutions  for  the  velocities  and  under  what 
conditions  the  more  involved  liquid-loaded  solid  formulism  must  be  used. 


It  was  found  that  there  are  certain  theoretical  similarities  between  the  surface 
waves  on  the  bulk  solid  and  on  the  plate  su-id,  which  were  pointed  out  in 


Theoretical  similarities  of  Rayleigh  and  Lamb  modes  of  (Q) 
vibration 


Pol«  of  the  infinite  plane-wave  reflection  coefficient  are  used  to  show  a correspondence  between  Rayleigh 
and  Lamb  modes  of  vibration  It  is  demonstrated  that  a Rayleigh  vibrational  mode  is  a special  type  of 
Lamb  mode  of  vibration.  Further,  it  is  shown  that  it  should  be  capected  that  one  vibrational  mode  for  a 
thick  plate  should  be  similar  to  the  theoretically  predicted  vibrational  mode  of  an  infinite  half  space,  a 
Rayleigh  mode.  Thus,  it  is  consistent  to  use  a thick  plate  as  an  approximation  to  an  infinite  half  space  and 
expect  results  predicted  by  Rayleigh-wave  analysis. 


n 


n 

Attention  was  then  focused  on  nonspecular  reflection,  since  it  was  found  that 
reflection  at  a Lamb  angle  or  at  the  Rayleigh  angle  is  accompanied  by  marked 
changes  in  the  profile  of  the  reflected  beam,  and  this  nonspecular  reflection 
in  turn  can  be  used  as  an  indicator  that  a plate  mode  has  been  set  up.  The 
theory  of  this  phenomena  was  given  in: 


Theory  of  Nonspecular  Reflection  Effects  for  an 
Ultrasonic  Beam  Incident  on  a Solid  Plate  © 

in  a Liquid 

Abstract- Pole*  ind  zeroei  of  the  infinite  plane  wave  amplitude 
reflection  coefficient  are  used  to  derive  a theoretical  prediction  of  the 
nonspecular  reflection  effects  which  have  been  observed  for  an  ultra- 
sonic beam  incident  on  an  isotropic  solid  plate  in  a liquid.  It  is  shown 
that  there  are  two  types  of  nonspecular  reflections,  and  they  can  be 
characterized  in  terms  of  a single  parameter  which  requires  a knowledge 
of  the  imaginary  part  of  a pole  of  the  infinite  plane  wave  reflection 
coefficient.  Theoretical  predictions  of  nonspecular  reflection  intensities 
are  presented.  Finally,  it  is  shown  that  the  reflection  characteristics 
for  high-frequency  beams  incident  on  thick  plates  are  the  tame  as  those 
expected  for  the  reflection  from  a single  liquid /solid  interface,  i.c.,  two 
infinite  half-spaces. 


i 


This  theory  enables  one  to  compare  reflection  from  hslf-spnces  ana  solid  plates. 
I was  found  that  non-specular  reflection  can  then  be  related  to  plate  thickness 
and  frequency,  as  was  shown  in: 


Non-specular  reflection  from  solid  plates 
and  half-spaces 


Recent  theoretical  developments  for  the  reflection  effects  seen  for  an  ultrasonic  bounded 
beam  reflected  from  a solid  surface  are  presented  and  a discussion  of  the  implications  of 
these  developments  is  related  to  the  case  of  reflection  from  infinite  half-spaces  and  solid 
plates. 


« 
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The  next  obvious  question  which  arose  was  concerned  with  the  relation  of 
reflection  to  transmission.  Although  this  relationship  does  exist,  a clearcut 
theoretical  analysis  is  rather  difficult,  and  thus  the  topic  was  not  actively 


followed  after  publication  of  the  following: 

Ultrasonic  bounded  beam  reflection  and  transmission  effects 
at  a liquid/solid-plate/liquid  interface 

Schlieren  techniques  are  used  to  show  that  bounded  beam  reflection  effects  occur  at  liquid/solid- 
plate/liquid  (L/SP/L)  interfaces  which  are  analogous  to' the  bounded  beam  reflection  effects  reported 
previously  at  liquid/solid  (L/S)  interfaces.  At  L/SP/L  interfaces,  nongeometric  effects  are  shown  to  be 
present  in  both  the  reflected  and  transmitted  beams  when  the  incident  angle  of  an  ultrasonic  beam 
corresponds  to  the  Lamb  angle.  In  addition,  similarities  between  wave  phenomena  at  L/SP/L  interfaces 
and  L/S  interfaces  arc  presented  which  suggest  that  the  description  of  bounded  beam  reflection  derived  by 
Bertoni  and  Tamir  for  a L/S  interface  can  be  qualitatively  applied  to  the  L/SP/L  interface  case.  This 
theory  is  shown  to  account  for  the  lateral  extent  to  which  a leaky  Lamb  surface  wave,  excited  by  mode 
conversion,  propagates. 


I 
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The  determination  of  plate  thickness  by  means  of  reflection  measurement  at 
onlique  angles  of  incidence  was  found  to  become  more  involved  whenever  the 
product  frequency  tine  thickness  became  larger  than  about  three  (in  units  of 
meter  tine  kHz) . Thus,  a study  was  undertaken  of  the  reflectivity  and  node 
structure  for  this  region,  contained: 


NON-SPEQULAR  ULTRASONIC  REFLECTION  FROM  THICK  SOLID 
PLATES. 

Non-specular  reflections  of  a bounded  ultrasonic  beam  from  thick  plates  in  water 
are  observed,  ft  is  shown  that  the  reflected  beam  profile  does  not  only  depend  on 
the  product  frequency  of  the  ultrasound  times  plate  thickness,  but  it  ultimately 
depends  on  the  imaginary  part  of  the  Lamb  pole  location.  Furthermore,  results 
indicate  that  a thick  plate  may  not  be  a valid  approximation  of  an  infinite  half- 
space, as  has  been  commonly  assumed. 


The  nonlinear  behavior  of  Loir.b  modes  on  plates  was  investigated  for  solid 
•nletes  in  order  to  detendne-  whether  or  not  a plate  can  indeed  support  a non- 
linear wave.  This  clearly  was  of  great  inportar.ee  since  the  actual  signals  used 
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in  the  Arctic  racy  contain  a number  of  frequencies,  and  it  was  not  at  all  clear 
whether  the  ice  plate  would  or  could  act  as  a frequency  filter.  The  basic  study 
of  the  behavior  of  the  plate  in  response  to  forced  nonlinear  oscillations  was 


reported  on  in : 
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NONLINEARITIES  IN  LAMB  WAVES  Qg 

The  equation*  of  motion  become  nonlinear  for  finite  amplitude  Lamb  modes  propagating 

In  Isotropic  plates.  Hence,  mode  Interactions  may  occur.  Two  such  Interactions  are 

the  Lamb  mode  three-phonon  Interaction  and  Lamb  mode  second  harmonic  generation. 

These  two  interactions  are  Investigated  experimentally. 

This  latter  contribution  discusses  basic  concepts  and  results,  and  it  together 
with  other  results  in  connection  with  nonlinear  phenomena  formed  the  basis  of  the 
renewal  proposed  to  Code  ^21. 

The  final  phases  of  the  Contract  work  were  concerned  with  the  ultimate  goal 
of  the  contrect:  the  study  of  the  ice  response  when  it  is  bounded  by  weter  on 
one  side  but  air  on  the  other.  Again,  the  investigation  of  the  magnitude  of  the 
reflection  coefficient  from  such  an  asymmetric  system  is  an  indicator  of  the  ice 
thickness.  Thus  the  formulism  for  this  system  had  to  be  developed  first  before 
the  already  available  experimental  data  (scaled-down  ice  pletes)  could  be  evalua- 
ted. The  theoretical  study  is  included  in  this  Report  in  its  entirely. 

Based  on  these  findincs  we  investigated  the  experimental  details  o:'  thickness 
determination  from,  the  reflection  of  bounded  width  ultrasonic  beams  from  different 
thicknesses  of  ice.  The  results  of  the  determination  of  the  node  structure  of  the 
ice  plates  are  now  being  compared  to  the  computation  model.  This  model  (making 
use  of  a computer  "Simplex"  nethod)  was  used  to  calculate  numerically  the  mode 
structure  of  the  system  air/ice  plate/sea  water.  A great  number  of  difficulties 
inherent  in  the  computer  nethod  made  the  conclusion  of  the  program  impossible 
within  the  estimated  tine.  However,  sore  of  the  preliminary  results  addressing 


the  eir/ice/water  system  have  been  published.  Concerning  the  group  and  phase 
velocity  of  the  surface  vave  created  in  the  process  of  critical  angle  reflec- 
tivity, a paper  vas  given: 


Sonic  phase  and  group  velocity  In  Ice  platea  Hotting  on 
water..  Peter  H.  Huang  and  Walter  G.  Mayer  (Physics  Depart- 
ment, Georgetown  University,  Washington,  DC  20057) 

. An  expression  for  the  phase  velocity  dispersion  of  vibrational 
inodes  for  an  Ice  plate  floating  on  water  Is  obtained  which  Is 
consistent  with  the  seemingly  different  expressions  contained 
In  the  reflection  coefficients  given  by  various  authors.  Ex- 
amples for  different  Ice  parameters  are  given,  and  the  disper- 
sion relations  for  the  group  velocities  Is  shown  to  be  more 
sensitive  to  variations  In  Ice  properties  than  is  the  phase  ve- 
locity. (Work  supported  by  the  Office  of  Naval  Research,  U.S. 
Navy.) 


More  recently  it  was  found  that  there  exist  severe  restrictions  on  the 
generation  of  the  leaky  Rayleigh  mode  when  a combination  of  water  and  ice  is 
used  as  the  boundary.  This  result  had  not  been  expected  and  was  found  on 


theoretical  grounds  very  recently: 


Investigation  of  the  conditions  for  the  existence  of  a leaky  Ray- 
leigh wave.  W.  G.  Mayer,  N.  G.  Brower,  and  D.  E.  Himbcrger 
(Department  of  Physics,  Georgetown  University,  Washington,  DC 
20057) 

The  Rayleigh  wave,  an  inhomogeneous  surface  wave,  exists  for  all 
isotropic  elastic  solids.  When  the  free  surface  of  the  solid  is  bounded 


by  a liquid,  a leaky  Rayleigh  wave,  or  inhomogeneous  damped  inter- 
face wave,  may  exist.  In  the  limit  that  the  density  of  the  liquid  goes 
to  zero,  the  leaky  Rayleigh  wave  tends  to  the  free  Rayleigh  wave. 
However,  the  leaky  Rayleigh  wave  does  not  exist  for  all  liquid/ 
isotropic  solid  systems.  A well-known  condition  for  the  existence  of 
the  leaky  Rayleigh  wave  is  that  the  velocity  of  sound  in  the  liquid 
must  be  less  than  that  of  the  shear  in  the  elastic  solid.  Upon  in- 
vestigation of  the  secular  equation  for  the  leaky  Rayleigh  wave 
velocity,  other  necessary  conditions  for  existence  become  apparent. 
The  conditions  arc  influenced  by  the  density  ratio  and  tbe  ratios  of 
the  various  velocities  in  tbe  liquid/solid  system.  [Work  supported 
by  ONR.) 


The  implications  of  this  last  finding  are  now  used  in  the  Simplex  computer 
method  to  arrive  at  meaningful  results  for  the  thickness  measurement.  However, 

I 

as  was  pointed  out  above,  this  new  analysis  is  not  complete  at  the  present  time. 
The  final  results  will  be  submitted  as  an  addendum  to  this  report  as  soon  as  the 


ir 


work  is  completed 


Volume  40 


ACUSTICA 

S.  HIRZEL  VERLAG  STUTTGART 


1978 


Plane  Wave  Reflection  from  a Plate  Immersed 
in  and  Floating  on  a Liquid 

by  P.  H.  Huang  and  W.  G.  Mayer 

Physic*  Department,  ( leorgetown  University,  Washington,  D.  20057 


Summary 

The  reflection  coefficient  for  a solid  plate  bounded  by  dissimilar  fluids  is  derived.  The  expression 
can  be  made  to  be  consistent  with  Brekhovskikh*  a. layered  system  reflection  coefficient  and  can 
be  reduced  to  the  reflection  coefficient  for  a two-lavered  system  derived  by  Schoch  and  to  the 
dispersion  equation  for  a liquid-solid-vacuum  system  given  by  Kwing.  Jardetaky  and  Press 
provided  one  takes  into  account  the  appropriate  forms  of  wave  propagation  and  coordinates 
of  the  systems. 

Reflexion  ebener  Welle n an  einer  beiilneillij  flueiiigkeilebegrenzlen  Piatu 
Zusammenfassung 

Ks  wird  der  P.eflexionsfaktor  einer  festen,  (lurch  unterschiedliche  Kluaaigkeiten  begrenr.ten 
Platte  hergelcitet.  Der  Ausdruck  kann  in  Ubercinstimmiing  mit  Brekhovskikhs  Keflex ionsfak tor 
fiir  ein  a-schichtiges  System  gebracht  und  auf  den  von  Schoch  abgeleiteten  Keflexionafaktor  fiir 
ein  xweischichtiges  System  zuriickgefuhrt  warden.  Beriicksichtigt  man  die  entsprechenden 
Formen  der  Wellenausbreitung  und  der  Systemkoordinaten,  so  kann  die  angegela-ne  Beziehung 
femet  auf  die  von  Kwing,  Jardetaky  und  Press  angegebene  l)is|M>rsionsglcichiing  fiir  ein  System 
Fluasigkeit-Festkdrper-Vakuum  zurOckgefiihrt  werden. 

Reflexion  d une  onrle  / Jane  /nr  une  /Jayne  immeryee  el  /lotion!  enlre  deux  liguides 
Sommaire 

On  calcule  le  coefficient  de  reflexion  d une  plaque  solide  illimitcc  sepunint  deux  fluides  diffc- 
rents.  I/ex  precision  trouvec  est  en  concordance'  avec  le  coefficient  de  reflexion  eleduit  du  systeme 
a a couches  de  Brekhovskikh  et  pent  etre  ra int'uc  egiilement  an  coefficent  trouvei  par  Schoch 
pour  un  svst^me  a deux  couches  ainsi  qua  Icquation  de  dispersion  pour  un  systeme  liquide- 
solide-vide  etablie  par  Kwing.  Jardetaky  et  Press,  pourvu  qu'on  donne  mix  propagations 
d'ondes  et  aux  coordonnees  dca  svstemes  lc-a  formes  convenables. 


1.  Introduction 

In  recent  years  a number  of  theoretical  and  ex- 
perimental studies  of  the  propagation  phenomena 
of  elastic  waves  in  a two-layered  or  a three-layered 
system  have  emerged  which  make  use  of  an  expres- 
sion for  the  reflection  coefficient  or  a dispersion 
equation  for  plate  mode  velocities.  A three-layered 
system  is  defined  by  a solid  plate  bounded  by  dis- 
similar fluids.  If  the  solid  plate  is  Ixiundcd  on  Ixith 
sides  hy  the  same  fluid,  it  is  calk'd  a two-layered 
system.  The  reflection  coefficient  is  referred  to  the 
amplitude  ratio  of  the  reflected  wave  with  respect 
to  the  incident  wave.  The  dispersion  equation  of  a 
layered  system  is  contained  in  the  denominator  of 
the  reflection  coefficient,  and  the  plate  mcxle  veloci- 
ties can  be  obtained  by  setting  the  denominator 
equal  to  zero.  Schoch  [ 1 1 obtained  an  expression  for 
the  reflection  coefficient  by  using  a particle-displace- 
ment formalism  and  applied  it  in  an  experimental 
study  [2]  of  a two-layered  system.  Ewing.  Jardetaky 
and  Press  [3]  gave  the  dispersion  equation  for  a 
floating  solid  plate  without  including  the  effect  of 


the  air  Ixiunding  the  plate  on  one  side.  Later  Brek- 
hovskikh  (4|  obtained  an  expression  for  the  reflec- 
tion coefficient  for  an  »- layered  system,  and  thus 
also  for  a three-layered  system,  by  using  a particle- 
potential  formalism. 

Unfortunately,  Brekhovskikh "s  expression  for  a 
three- layered  system  dex's  not  reduce  to  the  two- 
layered system  given  by  Sehcx-h.  The  differences 
occur  in  the  sign  of  the  imaginary  terms  and  there- 
fore the  hx'ation  of  a pole  of  the  reflection  coef- 
ficient in  the  wave-vector  plane.  Also.  Brekhov- 
skikh's  expression  for  the  dispersion  equation  for 
a three-layered  system  obtained  from  the  denomina- 
tor of  the  reflection  coefficient  dex's  not  reduce  to 
the  form  for  a liquid-solid-vacuum  system  given  by 
Ewing  et  al.  The  difference  appears  that  the  latter’s 
expression  does  not  contain  imaginary  terms  ex- 
plicitly while  the  other  formulas  do. 

It  is,  therefore,  the  object  of  this  article  to  resolve 
the  problem  of  discrepancies.  Extending  Schoch ’s 
treatment  to  derive  the  reflection  coefficient  for  a 
three-layered  system  Huang  (5)  had  taken  into  ac- 
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count  the  spatial  coordinates  with  respect  to  the 
selection  of  direction  of  wave  propagation  in  the 
layered  media.  It  is  shown  hen*  that  Brekhovslcikh's 
expression  of  reflection  coefficient  for  a three-lay- 
ered system  can  l>e  made  to  lie  consistent  with 
Huang's  expression  and  reduces  to  Schwh's  expres- 
sion for  a two-lavered  system,  and  that  the  expres- 
sion of  Ewing.  Jardetsky  and  Press  for  the  disper- 
sion equation  for  a liquid-solid-vacuum  system  can 
l>e  made  to  be  consistent  with  both  the  correspond- 
ing Rrekhovskikh's  and  Huang's  expression,  pro- 
vided one  takes  into  account  the  choice  of  the  co- 
ordinates system  and  appropriate  forms  of  wave 
propagation. 


2.  Derivation 

2.1.  The  derivation  from  Hrekhorslcilch's  expression  to 
Huang's  expression  and  the  reduction  to  Schoch's 
expression 

Assuming  a wave  propagation  in  the  form 

exp[  f i(fc  • * — m/)| , when'  x — x£hz£, 

and  using  the  spatial  coordinates  as  shown  in 
Fig.  la.  Schoch  [ 1 ] obtained  the  following  reflection 
coefficient  (disregarding  t he  phase  angle  tod  cos  fc\ ) 
for  a two-layered  system  with  plate  thickness  rf: 

* = (I) 

(.V  — iZ i sin  ysin  t ) ( Y — \Z\  cosy  cos  r)  -f  i.JZ\ 

( .V  — iZi  sin  y sin  f ) ( Y -f-  iZi  cos  y cos  t ) 

where 

.V  ~ cos*  (2  p)  cos  y sin  e -+-  4(c./c i)*  X 
X sin  x cos  x sin  (1  cos  p sin  y cos  t . 

>’  2 cos* (2 /f) sin y cose  + 4(c,/ci)*  x 
X sin  x cos  x sin  / 1 cos  p cos  y sin  r , 

J 2 cow*  (2  P)  sin  f cos  e + 4(c,/ci)*  x 
X sin  x cos  x sin  P cos  p sin  y cos  y , 

Zl  = Pi  ri  COSX/ppCi  cos  Oi 

with 

y 2 cudcosx/2ci, 
c h aid  cos  PI2ct. 

The  angles  ©i , x.  P are  connected  by  Snell's  law 
sin  0 1 sin  x sin  P 

(2) 

Cl  Cl  c» 

where  0i  is  the  incident  angle  in  medium  1,  a and  P 
are  the  angles  formed  by  the  respective  normals  to 
the  longitudinal  and  shear  wave  fronts  with  the 
positive  z-axis,  Ci  is  the  wave  velocity  in  the  liquid. 


2 


Fig.  1.  Spatial  coordinates  with  respect  to  the  incident  and 
the  reflected  plane  waves  used  by  (a)  Schoch  for  a two- 
media  system,  (b)  Rrekhovakikh  for  a three-media  system, 
and  (c)  Huang  for  a three-media  system. 


ci  and  c,  are  the  respective  longitudinal  and  shear 
wave  velocities  in  the  solid  plate. 

Considering  the  same  waveform  but  a different 
coordinate  system  with  respect  to  the  incident  and 
the  reflected  waves  as  shown  in  Fig.  1 b,  Brekhov- 
skikh  (4]  obtained  the  following  reflection  for  a 
three-layered  system,  based  on  his  expression  for 
an  n- layered  system 

(3) 

= M (Z,  - Zs)  + i[(**  - M*)Zt  - Zil 
M (Zi  + Z,)  + i[(N*  - M*)Z,  + Z,]  ’ 

where 

Z2cos*(2y2)  Z«sin*(2y2) 

Z,sinP  + Zi  sin  Q ’ 

_ Z2cos*(2y2)cot  P Z2(sin*(2y2)cotQ 

Mm  * + IT  * 

Z i = pici/cos  ®i , Zf  = piCt/oos  02 
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Z>  o2 c a/csss  O i . Z-i,  ii  i 62/cos  y> , 

/’  in  d com  Oijc 2 , Q tad  cos  y2/02 . 

The  corresponding  Snell's  law  assradated  with  eq.  (3) 
to  eq.  (2)  is  as  follows: 

sin  0i  sin  02  sin  y2 


Extending  Schoch’s  work.  Husnu  |5|  derived  the 
following  reflection  coefficient  for  a three -layered 
system  as  shown  in  Fig.  1 c by  taking  into  account 
the  positive  ^-direction  to  lie  either  upward  or 
downward  and  the  wave  propagation  to  Is-  either 
of  the  form 


exp[4- i(fc  • * — c»f)|  or 
exp[ — i(k  ■ x — oil) | , 


when-  x — xi 


By  adding  and  subtracting  some  identical  terms, 
the  denominator  of  eq.  (7)  liecomes 

2|(Zi/Z2)2(.V*  - .1/2)  + i.W(Z,/Z22)  x 
x(Z,+Zs)-rZsZ,/Z2*|  = 

= 2|<Z,/Z2)*JV*  - (Z,/Z2)*J/*  - 
- (Z,/Z2)2  MS  + (Z,/Z2 )*  MS  + 

4-  i(Z,/Z2)*  M 4-  ifZjZi/Zj1)  M + 
f i (Z j /Z2)2.V  — i (Z3Z1/Z22)  S 4- 
-f  ZsZ,/Z-.2|-2i(Z,/Z2)*^  + 

+ 2 i(ZsZi/Z22)  .V  = 

- 2(Z,/Z2)[AT  + M - i(Z,/Z,)]  x 
x (Zi/Z->)  (AT  — M fi)- 
-2i(Z,/Z,)X[(Z,/Z2)-(Z,/Z2)|.  (8) 

Similarly,  the  numerator  of  eq.  (7)  can  la-  written 
as 

2[(Z,/Z2)*(JV*  - .1/2)  + i.W(Z,/Z22)  x 
x (Z2  — Zj ) — Z3Zi/Z22)  = 


(2X4-2  iZs  sin  y sin  r)  ( )‘  4-  iZi  cos  y cos 
(2X4-2  i Z3  sin  y sin  f)  ( )'  ^ iZi  cos  y cos 


r)  ± i J(Zt  + Zs) 
F)  "f  iJ(%\  — Z3) 


Z3  03  rs  cos  3.) iip  ri  cos  03  . 

03  is  the  angle  lietween  the  normal  to  the  trans- 
mitter! wave  front  in  medium  3 mid  the  positive 
:-axis;  X.  }'.  J.  Zi . y and  f are  defined  as  in  eq.  (1 ). 
Choosing  the  positive  i-direetion  to  la-  upward,  the 
upper  sign  in  eq.  (5)  should  la-  used  if  the  wave 
propagation  in  the  media  is  descrilad  by 

exp|  -f  i (Ir  • x — ml)] , 

and  the  lower  sign  if  the  wave  is  given  by 
exp[—  i (Ar  - jr  — nif)| . 

Reversing  the  direction  of  the  poaitivc  r-axis.  the 
upper  sign  should  be  used  in  eq.  (5)  if  the  wave  is 
given  by  the  negative  exponential,  and  the  lower 
sign  for  the  wave  descrilaal  by  the  positive  ex- 
ponential. 

Following  this  rule.  eq.  (3)  can  la-  rewritten  as 
shown  in  eq.  (fi)  under  the  transformation  of  the 
physical  system  pictured  in  Fig.  I b to  that  in 
Fig.  lc, 

M (Z3  — Z,)  — i[(iV*  — M*)Z\  — Zt] 

M (Zj  -|-  Z\)  — i[(JV*  — ;W*)Zi  -f  Zs|  ( ' 

By  multiplying  laith  the  denominator  and  the  nu- 
merator by  a common  factor  2i(Z|/Z22),  eq.  (3)  can 
la>  written  as 

V = (7) 

2(Z,/Z22)  [i  A/(Z,  - Z,)  + (A'2  - Jt/*)Z,  - Z,1 
“ 2 (Z1/Z22)  [i  M (Z,  + Zi)  -+-  (2V*  — M*)Z\  + Z*)  ’ 


= 2(Z,/Z2)[X4-  M - i(Zs/Z,)]  x 
x (Z,/Z2)(X  - M -i)  + 

4-  2i(Zi/Z2) .V|(Zj/Z2)  4-  (Zs/Z2)|.  (9) 

Now 

(Z,/Z2)|X  4-  M -i(Z,/Z,))  = 

[cos*(2y2)  cos  ( /*/2)/sin  ( /*/2)]  4- 
4-  |(62/r2)  (cos  Oil  cos  y2)  x 
x sin*(2y2)  cos(Q/2)/sin(Q/2)]  — 

— i 03  C3  cos  02/«2  r2  cos  03 , (10) 

(Z,/Z2)(.V-  M ± i)  = 

|cos'2(2  y2)  sin  ( />/2)/cos(/,/2)j  4- 

+ [(6*/e*)  («**  02/cos  y2)  x 

x sin*(2y2) sin(Q/2)/cos(Q/2)|  ± 

rfc  ipi  fi  cos  02/p2e2  cos  0i , (11) 

and 

2i(Z,/Z2)AT  = 

= i{[cos2(2y2)/sin(/,/2)cos(/,/2)]  + 

4-  (02/f2)  (cos  02/cos  y2)  x 
x sin2(2y2)/sin(^/2)  cos(Q/2)]} . (12) 

By  making  a change  of  the  following  corresponding 

notations  used  in  eq.  (3)  and  (5), 

P/2  -*■  y . Ql%  -*■  e . 

c*  -*•  Ci . 02  -+  « , £>*  -*  <?p  • 
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E<|s.  (10).  (11)  and  (12)  become,  respectively, 

(Zi/Z2)[(iV  + V - i(Zj/Zi)|  = 

— (cos2  (2/1)  cos  y/sin  y)  f4(r„/f|)'2  X 
X sin  a cos  * sin  p cos  0(cos  r/ sin  e ) — 

— iZs  = (X/sin  y sin  e)  — i Z3,  (10') 

(Zi/Z-i)  (N  - -»/  ± i)  = 

= | cos2 (2//)  sin  y/cos  y|  -(-  4(r„/e i)2  X 
X sin  a cos  a sin  p cos  p (sin  e/cos  <■•)  4_ 

±iZ,  = (r  /cos  y cos  f ) i iZt  (11’) 

and 

2i(Zi/Z-.)A’  = i {[cos2  (2  P)l»\n  y cos  y | + 

-f  4 (r,/ci)2  sin  a cos  a sin  p cos  pi  sin  t cos  t [ 

= iJ/siny  cosy  sine  cos  r,  (12') 

where  Snell’s  law.  eq.  (4).  has  Iteen  used. 

Substituting  cqs.  (10').  (11')  and  (12')  into  eq.  (R). 

one  obtains  for  the  denominator  of  eq.  (7) 


vacuum  system  pictured  in  Fig.  2 where  the  wave 
propagation  is  of  the  form  expl  f i(«il  — *•*)! 

P(2Q  -(-  d cosh  mi  //  cosh  vi  H)  -t- 

— y A sinh  Vi  H sinh  v\  H ~ 0.  (15) 

Eq,  (15)  can  la*  rewritten  as  follows: 

2(PQ  f Pd  cosh  ri  //  cosh  vi' H)  4 
-f - d(  — P cosh  vi  II  cosh  vi  H t 
4-  Q sinh  vi  H sinh  »i'H)  =0,  (10) 

where 

/>  = 4_  jfct)2  cosh  »i  II  sinh  v\  H — 

— 4 ri  v\  k*  sinh  »i  H cosh  v\  H , 

Q (i<i'2  t k-)'1  sinh  »*i  //  cosh  v\  H — 

— 4 i*i  vi'k * cosh  »'i  II  sinh  v\  H , 

d Q2*2i(vi'i  - **)  (1-22  - F2)  vi/pi  pr  1*2 


2|(Z,/Z2)2(.V2  - .l/2)  4 i J/(Z,/Z2-)(Z,  4 Zs)  -+  Z,Z,/Z22|  = 

(2  X — 2 iZs  sin  y sin  r)  ( Y 4 iZi  cos  y cos  e ) — i J (Z\  — Z3) 
sin  y cos  y sin  f cos  f 

Substituting  cqs.  (10’).  (11’)  and  (12')  into  eq.  (0).  one  obtains  for  the  numerator  of  eq.  (7) 

2l(Z,/Z2)2(N*  - m 4-  i.V(Z,/Z**)(Z3-Z,)  - Z3Z1/Z,2)  = 

(2  X — 2 tZa  sin  y sin  e)  ( >'  — iZt  cos  y cos  f)  4 \J  (Z,  4 Z3) 
sin  y cos  y sin  f cos  f 


(13) 


(M) 


Finally,  by  dividing  eq.  (14)  by  eq.  (13).  one  thus 
obtains  exactly  the  same  expression  as  eq.  (5)  in 
accordance  with  the  rules  concerning  the  sign  of  the 
maginarv  terms  as  discussed  previously. 

Thus.  Brekhovskikh’s  expression  eq.  (3),  is  con- 
sistent with  Huang’s  expression,  eq.  (5),  provided 
one  takes  into  account  the  coordinates  of  the  sys- 
tem with  respect  to  the  incident  and  the  reflect'd 
waves. 

Next,  as  a special  case  when  the  solid  plate  is 
bounded  on  laith  sides  by  the  same  liquid,  i.c..  a 
two-layered  system  as  shown  in  Fig.  la.  the  nor- 
malized impedances  of  medium  1 and  medium  3 
are  the  same,  namely.  Z3  = Z\.  It  is  now  evident 
from  the  derivation  shown  alsive  that  lnith  eqs.  (3) 
ami  (5)  reduce  to  an  expression  which  is  identical 
to  Schoch’s  expression,  eq.  (1). 

2.2.  The  reduction  of  Ewing'.  Jardetxky  and  Pres* 
expreexion  to  Hrekhovekikh'e  and  Huang  * ex- 
prtssion 

Eq.  (15)  gives  Ewing.  Jardetaky  and  Press’  ex- 
pression [3|  of  dispersion  equation  for  a liquid-solid- 


wit  h 

v,  (t*-*2.,)''*- 

.-2  5(F2  — *=.)'  *. 

The  corresponding  notations  used  in  eqs.  (IB)  and 
(5),  respectively,  are  as  follows: 

p2  C>1  • PI  **  ?P  ' *1  ■**  C|  . 

*2-**'Cl.  /fl**-Cg,  //»*rf/2. 


O tax’! 


Fig.  2.  Spatial  <*oordinates  used  by  Ewing.  Jardetaky  and 
Prnw  for  a liqiiid-aolid-vacuum  ayntem. 
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n •»  — ioi con */f| , n'  ••  — i <i> cos /#/«•„ , 

»•*  **  — I frj  COS  ®|/C|  . 

It  is  noted  that  Snell's  law  is  associated  with 
eq.  (Hi).  and  its  connection  with  angles  x.  fi  and 
(t\  in  eq.  (5)  viz., 

k — k%l  sin  x <■,(,  sin  fi  sin  (i\ . 

has  Isssi  used  to  ohtain  the  corresponding  terms  i'i . 
r i'  and  v>.  Henee.  in  eq.  (Iti).  one  obtains 

eosh  Vi  II  |exp(  icirf  cos  x/2f|)  4 

4 exp(  (-  iinrf (sis x/2ri)|'2 

cos  (did  cos  x/2r|)  = 

cos y . (17) 

cosh  i'i’//  |exp(-  ioirf cos  fil'2r,)  • 

- cxp(  f-  i*n d cos  [i  2r,)|  2 
cos((i»d  c< is  fi  ‘2 r„) 

cos  f , (IK) 

sinh  ri //  |exp(  i tod  cos  x 2 <-i ) 

exp(  ■ icirfcos  x 2f|)|  2 
isiny.  (lit) 

siidi  v\  H |exp(  — i tod  cos  (i!‘2c,) 

cxp(  • Mod  <sis  fi  2r„)|  2 

— isinr.  (20) 

( r »’*+**)  (2 **  - k", )-’  Jrj,  co*-(2/() . 

(21) 

— 4 ri  V| ' i-'-  - 4(  imcosxri)  ■ 

x ( — in>  i'os  fi  r,)  (<»  sin  x,V|)  («i  sin  fi  r,) 

4 k*t  (r.  'C|)z  sin  x cos  x sin  (i  cos  fi , (22) 

A |«l(-  i«H  COS  X C|)  «,,(  i(K  COS  0)  C))|  ■ 

X (»'l/<,s)*(^l ) (r**/*-!  )*  = 
k\ , pi  f|  cos  x/ppr | cos  0, 

= KA  (23) 

and 

/’  - — i/rj,  |eos2(2//)  cos  y ns  r ■ 

4 4(r,/ci)*  sin  x cos  x sin  fi cos  fi  • 

X sin  y is**  f|  = 

= -'*i,X.  (24) 

V " i/tj,  |(sis2(2//)  sin  yens  f * 

4 4(r„/r i)1  sin  x cos  x sin  fi  isis  (i  x 
X rsm  y sin  f | = 

— »*;.  r.  (25) 


where  A'  and  >'  arc  as  defined  in  eq.  (5). 

Now 

pg  = (-i*S1JD(-i*j1  r)~-k$,xr, 

/'A cosh  i'i  //  cosh  ri'//  = 

= (-  >*i,  X ) (tj,  Z\ ) cos  y cos  e = 

— it^,  A'Xi  cos  y cos  r (27) 

and 

— /'cosh  vj  II  cosh  v\  H 4 
4 V sinh  i’i  //sinh  v\  H = 
i(-Jl|<os2(2//)  cosy  sin  f • 4(f„/ci)2  x 
x sin  x (s>s  x sin  fi  cos  fi  sin  y cos  f|  X 
X cos  y cos  f 4-  | csis2  (2  fi)  sin  y cos  f 4 
• 4(f,/ri)s  sin  x cos  x sin  fi  cos  fi  x 

x (sis  y sin  f|  sin  y sin  f = 

ilrj,  | (sis'- {‘2  fi)  sin  f cos  f (cos2 y 4 sin2y)  • 

4 4(r„/ri)2sin  x (sis  x sin //cos//  x 
X sin  y (sis  y (cos2  f 4 sin2  f)  ! — 
ii-Jl|cos2(2//)sin  f cos f 4 4(r»/ri)2  x 
X sin  a (sis  x sin  fi  is  is  fi  sin  y cos  y | — 

= i k},J.  (2K) 

where  ./  is  as  dctiiiisl  hi  eq.  (5). 

Kulstitntini!  eqs.  (2ii).  (27)  and  (2H)  into  is|.  (Iti), 
one  obtains 

2 ( k't X )■  i i-J,  XZ,  (SIS  y (Sis  F) 

i*$,./2,  =0.  (2ft) 

Sims'  /•)),  4 tt,  thus  e(|.  (2ft)  Isssnncs 

2 -V  ( 1'  \Z\  cos  y cos  r)  — iJZ\  = 0 . (30) 

Next,  (sinsidcr  (S|.  (5)  for  the  case  where  nnsliuin  3 
is  a vacuum,  i.e..  0.  It  follows  from  (sp  (5) 

that  the  reflection  coefficient  for  the  case  where 
the  wave  is  of  the  form  exp|  t i (fol  — k ■ *)|  and 
the  s|iatial  (ssirdinates  as  shown  in  Kig.  2 can  lie 
written  as: 

2 A'  ( > i^i  (sis  y (sis  f)  4 \JZ , 

2 A”  ( Y 4 iZi  cos  y cos  f)  — iJZi 

The  dispersion  (spiat ion  for  the  liquid-solid-vacuum 
system  is  therefore  given  by  the  vanishing  of  the 
denominator  of  eq.  (31 ).  leading  to  exactly  the  same 
expression  bs  eq.  (30). 

It  has  Iss-n  shown  that  Kwing.  .lardetskv  and 
Press'  expression.  eq.  (I«).  is  indissl  consistent  with 
eq.  (31).  provided  one  takes  the  (ssirdinates  of  the 
system  and  the  forms  of  wave  propagation  into 
account. 


I 


228 


P.  H.  HI  ANG  and  W.C.  MAYER:  REFLECTION  FROM  \V  IMMERSED  PLATE 


u I’STICA 
v«i. ««  di;m 


Because  Brekhovskikli's  expression.  eq.  (3),  is 
consistent  with  Huang's  expression,  eq.  (.5),  there- 
fore Kwing.  Jardetsky,  and  Press'  expression, 
eq.  (lti),  is  also  in  agreement  with  Brekhovskikli's 
expression  for  a liquid-solid- vacuum  system. 
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3.  ('oncltisioil 

In  this  article,  we  have  reviewed  various  expres- 
sions of  reflection  coefficient  and  dispersion  equa- 
tion given  by  Schocli  [||,  Brekhovskikh  [4|.  Kwing 
et  al.  |3|  and  Huang  |5|.  The  result  is  established 
that  it  is  possible  to  make  these  expressions  consist- 
ent with  one  another  provided  one  takes  into 
account  the  coordinates  of  the  physical  systems  and 
the  appropriate  forms  of  wave  propagation. 
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